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This  report  studies  radar  backseat  ter ing  from  a  vegetated 
terrain.  Tne  vegetation  is  modelled  by  spherical  water  droplets 
which  car,  be  treated  «b  discrete  sea 1 1 ere r s  .  The  vegetation  is 
assumed  sufficiently  lossy  so  that  the  underlying  ground  is  not 
noticeable .  The  method  of  ?o i dy  is  used  to  evaluate  the  mean 
fieJJ  m  the  vegetation  when  the  wavelength  of  the  incident  radi 
*a.a vH  is  large  compared  to  the  droplet  size .  Once  an  equivalent 
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dielectric  constant  for  the  vegetation  is  obtained  by  the  Foldy 
technique,  single  scattering  is  employed  to  evaluate  the  back- 
scattering  cross  section.  The  resulting  expression  is  found  to 
compare  favorably  with  experimental  data.  In  addition  to  the 
three  dimensional  work,  a  one  dimensional  problem  is  analyzed. 

This  analysis  is  compared  with  the  Foldy  approximation  in  the  high 
density  limit. 
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I .  INTRODUCTION 


This  report  studies  microwave  backscatter ing  from  vegeta¬ 
ted  terrain.  Discrete  random  media  theory  is  employed  in  an 
effort  to  relate  the  backscattering  cross  section  to  charac¬ 
teristic  objects  in  the  vegetation.  The  resulting  theory  is 
then  compared  with  results  obtained  by  more  empirical  methods. 

The  study  is  motivated  by  the  need  to  relate  radar  return 
to  the  characteristics  of  the  scattering  objects.  This  situ¬ 
ation  is  particularly  complex  in  the  case  of  vegetated  terrain 
which  consists  of  an  ensemble  of  many  highly  irregular  objects 
placed  xn  a  more  or  leas  random  fashion.  Using  the  scattered 
return  fro®  a  vegetated  terrain  one  would  like  to  obtain  such 
information  as  height,  density  and  moisture  content  of  vege¬ 
tated  regions  in  addition  to  obtaining  the  characteristic 
shapes  of  scattarers,  so  that  the  vegetation  can  be  classi¬ 
fied.  This  type  of  information  would  be  important  for 
military  analysis  of  the  terrain.  In  addition,  if  the  terrain 


could  be  determined  adequately  from  an  eiectr 
of  view,  the  bacxscattsrer  information  could 
diet  returns  at  other  angles  and  frequencies 
This  would  be  of  great  aid  to  the  development 
data  bases  for  radar  simulations. 

The  above  applications  '  serv-  d  as  a 
the  development  of  slactrosiag.vt.c  models  for 
terrain .  These  models  have  bean  cons  ;.t  acted 
vegetation  with  a  s lab  of  random  medium 
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char.scte:  istics  a r®  related  to  the  physical  quantities  in  tne 
medium.  The  models  can  oe  divided  into  two  catagones: 
continuous  and  discrete.  In  the  continuous  case,  the  vegeta¬ 
tion  is  modeled  by  assuming  its  permittivity  c { r ) ,  is  a  random 
process  whose  moments  are  known.  The  average  b&ckscattering 
cross-section  is  then  ca  culated  fx'om  a  knowledge  of  the  sta¬ 
tistics  of  c.  (r)  .  Usually  just  the  mean  and  correlation  of  the 
permittivity  are  required.  Following  thus,  some  quasi  quanti¬ 
tative  techniques  are  used  to  relate  the  medium's  statistics 

to  the  actual  vegetation  under  cons ideration .  General  devel- 

( 1 ) 

optnent  of  these  techniques  are  attributed  to  Keller  and 

( 2 ) 

Tatarskii  and  Gertsenshteir.  ,  however,  particular  applica¬ 


tions  to  vegetated  terrain  have  been  made  by  Lang  , 

(4)  ,  „  (5) 

Hevenor  and  r ung 

When  modeling  is  done  by  the  discrete  technique,  on  the 
other  hand,  the  individual  objects  ir.  the  medium,  such  as 
leaves,  are  identified  by  their  determini s tic  cross  section, 
and  then  each  object  is  placed  randomly  and  its  position  is 
given  by  means  of  a  probability  density  function.  The  average 
backscattered  cross  section  is  computed  m  such  a  **ay  that  all 
multiple  interaction  between  objects  is  accounted  for.  The 


systematic  developmen 


these  techniques  has  been  made  by 
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Fold-/  ,  Lax  '  and  Twer sky 
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,n  the  continuous  case,  since  this  has  already  been  included 


urnmnuffwimw  f9»S»S^sS*SK<««aWW!WSP8!^SWWiJP*(PHi^^Wi 

. . .  iWV1***  i V 


As  mentioned  earlier,  *j«  will  employ  discrete  scattering 
techniques  to  model  the  vegetation.  We  assume  a  lossy  half 
space  of  vegetation  is  present,  i,e.„  one  in  which  the  ground 
is  not  visible.  We  than  consider  the  vegetation  to  be  made  up 
of  small  water  droplets  which  are  randomly  distributed .  In 
Section  II,  the  average  backscattering  cross  section  is  found 
by  assuming  the  electromagnetic  field  can  be  treated  as  a 
scalar  field.  The  Foldy  technique  is  employed  to  calculate 
the  field.  This  technique  requires  that  the  scatterers  be  in 
the  Rayleigh  regime,  i  .  ,  small  compared  to  the  wavelength. 

The  use  of  the  scalar  assumption  introduces  difficulties  in 
explicitly  evaluating  the  cross  section.  These  difficulties 
are  removed  in  Section  III  where  the  full  electromagnetic 
problem  is  treated  for  horizontal  polarization .  In  Section 
IV  numerical  results  are  presented  and  compared  with  data 
from  <-vther  ,  Certain  rest.  v.  lions  are  placed  on  the 

use  of  Foldy' s  technique  in  Section  V.  There  a  one  dimension- 
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is  analyzed 

exactly  and  examined 

in  the  thermodynamic 

limit.  A 

comparison 
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method , 
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II. 


ISOTROPIC  SCATTERING  BY  A  HALF-SPATE  OF  RANDOMLY 


DISTRIBUTED  SPHERES 


In  this  chapter  the  problem  of  multiple  scattering  of 
scalar  v/aves  by  a  half-space  of  randomly  distributed  dielec¬ 
tric  spheres  is  considered.  The  spheres  are  assumed  to  ha  of 
the  same  site  and  their  radius  is  taken  much  smaller  than  the 
wavelength  of  the  incident  wave.  In  addition,  we  also  assume 
that  the  spheres  are  independent  and  uniformly  distributed. 

An  approximate  expression  for  the  mean  field  is  derived 
by  employing  Eddy's  method^.  This  approximate  mean  field 
i*  then  used  along  with  the  Born  approximation  to  calculate 
tne  backscattered  field.  The  scattering  coefficient  is  then 
obtained  from  tne  backscattered  field. 

2-1  SINGLE  SCATTERING 

Consider  a  half-space  (2>0)  of  randomly  distributed 
spheres  (see  Fig.  2-1!.  A  scalar  wave  .^(r)  *  exp  (-jk^-r) 
with  harmonic  time  dependence  exp  { jwt.)  is  incident  upon  the 
scatterers .  The  angle  of  mc’  ience  is  r  .  In  free  space  „.,(r) 
satisfies  the  wave  equation. 

Y 3  if  -  <  xr )  *  k i  (  r )  -  0  ( 2-1 ) 

0  —  ')  u  — 

where  *  w»  uncn  and  and  c.  are  the  free  space  permaabil- 
tty  and  permit! vity  respectively. 

The  wave  scattered  by  a  single  sphere  which  is  small  com¬ 
pared  ?  •  •  wace length  is  given  by 


(2-2) 
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Here 


gQ(r)  «  exp  (-jk0  ~  E 


:  2-2*) 


is  the  free  space  scalar  Green’s  function  which  satisfies  the 
wave  eq uation : 
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The  scattering  p.cperties  of  the  spheres  are  characterized  by 
the  relationship 


A  -  v  v  1  i  r  i 

n  -n 


(2-4) 


which  makes  the  strength  of  the  scattered  wave  from  the  n— 

f  r*  } 

sphere  proportional  to  the  externa,  field  ,  r  >  acting  on 

— n 

it.  The  proportionality  constant  •  is  cal .ed  the  scattering 
coefficient  and  is  the  same  for  all  spheres.  The  simple 
spherical  wave  behavior  of  the  scattered  field  results  from 
the  assumption  that  the  fields  are  approximately  constant  in 
the  vicinity  ,>f  the  sphere.  As  a  result,  each  sphere  can  be 
replaced  by  a  point  source  located  at  its  origin  (Rayleigh  as- 
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-2  MEAN  FIELD 


Mow  consider  scattering  by  N  spheres  which  are  all  smal. 
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particular  scattercr  is  located  in  the  volume  V  is  independent 
of  the  locations  of  the  other  scattorers,  thus  v*  can  write 


P  ^  ^2, '  ii 2  y  “  *  *  ^  m  p  ( )  P  ( 1 2 )  •  •  •  P  ■  / 


where  p  ( r  ) 
r  — n 


1 

?* 


(2-5) 

(2-6) 


The  fundamental  equations  of  multiple  scattering  will 
rsov  he  formulated.  Consider  a  particular  configuration  of  the 
spheres.  Th©  wave  function  at  the  point  r  is  given  by 
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The  equations  (2-7)  and  (2-8)  represent  the  fundamental  equa¬ 
tions  of  multiple  scattering. 

Lvmg  the  problem  would  then  ccn- 
.nul taneous  linear  algebraic  Eqs , 
and  substituting  these  in  Eq .  (2-7), 

function  of  the  positions  and  scattering 
parameters  of  the  spheres.  Taking  the  mean  value  of  this 
quantity  would  then  give  us  the  desired  results. 


The  direct 

meth.cd 

c  t 

SI  St 

of.  solving 

th  e  s  e 

*  ^  4 

O  -J  *. 

(2-8! 

for  tne  ' 

n ) 

!  <  r  ? 

and 

thus 

x  *  i.  n  ^  *•'  \  x*  / 

as  a 

t'n  fortunately,  it  Joes  not  seem  possible  to  carry  out  this  pro¬ 
cedure  because  of  its  complexity  when  N  is  large  and  it  in 
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-n~l 1 


•  ;  r)  '4* (r  )dr,  , 
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The  quantity  <*(r  'r  )  >  represents  the  external  field 

— n  — n 

th 

acting  on  the  n  sphere  averaged  over  all  possible  configu¬ 
rations  of  all  the  other  spheres.  For  large  M  we  approximate 

th 

the  external  field  on  the  n  '  sphere  by  the  mean  field  which 

c  h 

would  exist  at  the  position  of  the  n  4  sphere  when  the  sphere 
is  not  present.  Thus  we  can  write 


< *  ( r  r  )  >  ■-  <*(r  )> 

-•n  — n  — to, 


(2-14) 


iy  using  Eq» .  (3-14)  ,  *1-13)  ,  Eq.  (.1-12}  .becomes 
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The  quantity  <:•*>  (  r  Vpy  r-r^  !  N  is  the  same  for  all  the 

spheres.  Substituting  Eg.  (2-15)  into  Eq.  (2-9;  we  obtain 
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But  from  Eq.  <2~6i  we  have  o-r  ’  “l.V  thus 
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To  see  t h«  physical  significance  of  the  Eq.  (2-18),  let 
u a  apply  the  operator,  7 1  4  ic*  to  both  sides  of  Eq .  (2-18) 


7s<v(r)>  4  k*<^(r)  >  ««  (7  *-*■**)  jQ  (r) 


Y.Jdr^lT^k'r^r-^)  <*(£„)> 


(2-19) 


By  using  Eqs .  (2-1)  and  (2-3),  Eq.  (2-19)  becomes 


’  *  <  •£  ( r )  >  4  k *  <  *  ( r )  >  -  0 
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where 


(2-20) 


k  *  »  k,t  f  jv 
eq  0 


(2-21) 


We  see  that  <v(j)>  in  the  region  z>0  satisfies  the  wave  equa¬ 
tion  in  a  "continuous  medium"  m  which  the  wav©  number  depends 
upon  the  scattering  coefficient  and  density  distribution  of 
the  spheres  . 

The  problem  of  finding  the  average  value  of  til*  wave 
function  has  beer,  essentially  reduced  to  solving  a  boundary 
value  problem  for  the  wave  equation  arid  determining  the  value 
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Carrying  out  this  calcu 


on,  we  find  that  the  mean 


wave  in  the  equivalent  or  effective  medium  is 
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x  eq  x  st¬ 

and 

k  *  Jcnsind  (2~26) 
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2-3  CORRELATION  FUNCTION  OF  THE  SCATTERED  FIELD 

We  now  obtain  an  integral  representation  for  the  total 
scattered  intensity.  The  second  term  of  the  right  hand  aide 
of  Eq.  (2-7)  gives  the  total  scattered  wave. 

M  ,  . 

*  <r)  **  :  vV'n;  Ir  }  g(r~r  )  (2-27) 
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Since  the  scattered  field  radiates  into  the  equivalent 
medium,  k  ,  w«  have  used  the  Green's  function  g(r>  instead 
of  g^lr)  whore 

-  jk  . ri 
J  eq  — 

g  (  r )  -  — — — —  12-23} 
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To  obtain  the  mean  value  of  : *  (r) 1 J  we  follow  an  analogous 

A  “* 

procedure  to  that  used  for  <  *(r) > ,  but  the  analysis  now  be- 
coi»<# a  somewhat  more  involved.  To  begin  we  first  multiply  the 
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We  next  assume  th* t.  the  external  field  at  the  n  sphere 

( 6 ) 

can  be  replaced  by  the  mean  field;  thus  we  have 
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If  we  let  h(r  .r  )  "9(£i-r)  g*  (r  5-*r)<*(r  )><y*(r  )>,Eq.  (2-31) 
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ia  written  as 
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Using  Eqs .  (2-5)  and  (2-6)  we  obtain  for  th*  bracketed  term 

at  the  right  hand  side 
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The  first  integral  consists  of  the  (N  -N)  t tarns  where  nj*m, 
while  the  second  integral  consists  of  M  terms  where  n*»m. 


Since 
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where  for  large  N  and  ^ 
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Substituting  Eq .  {2-  36)  into  Eq .  (2-  32)  we  have 
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2-4  EVALUATION  OF  THE  TRANSVERSE  SPECTRAL  DENSITY  AND  BACK- 


SCATTERING  COEFFICIENT 

In  this  section  we  evaluate  the  transverse  power  spec¬ 
tral  density  of  the  fluctuating  portion  of  the  scattered 
field.  From  this  spectral  density  a  direct  calculation 
yields  the  backscattering  coefficient  as  we  show  in  Appendix 
B. 

We  stati  by  defining  the  fluctuating  portion  of  the 
scattered  field  by 

■  v  (r)  -  <■-»  At)  >  (2-40) 

where  we  note  that  ■  ;>(r)  >«*0.  The  spectral  density  can  be  ob¬ 
tained  by  Fourier  transforming  the  correlation  function  of 
Vf(r),  We  proceed  by  calculating  the  correlation  of  ^ { r) : 
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Equation  (2-41 i  shows  that  at  r^*r7  the  fluctuating  in¬ 
tensity  is  the  total  intensity  minus  the  coherent  intensity. 
By  using  Eqa .  (2-39)  and  (2-10)  in  Eq.  (2-41),  we  find 

’  *  ~1 5  *'*  !  --2  !  >  "  g  (  *  )  g*  ( r^-r  ’ )!  •  ^  (£  ‘ )  >  i  2d£’  (2-425 

Mow  proceed  by  taking  the  transverse  Fourier  transform  of  Eq. 
(2~42)  with  respect  to  r  and  r  where 

C  ^  v  q 
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i-1*  2 


with  a  and  a,,  being  unit  vectors  in  the  x  and  y  directions 
-~x  — y  3 

respectively . 

W’<a  have 
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Here  we  have  defined  p^  by 
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and  k  *k  a  -*-k  a  , 

-t  x-~y  y  y 


We  now  *at  z, «z„-0  since  it  is  the  spectral  density  at 

X  «u 

the  interface  that  is  related  to  the  backscattering  coeffici¬ 
ent.  To  simplify  Eq,  (2-42)  represent  tube  Green's  functions 
by  their  transverse  Fourier  representation,  i.e. , 


*3 {<* :  z'*±t  z] 
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(2-45 


with  <  -*  x*  -  ■  k  * ,  and  use  the  explicit  expression  given  in 
z  eq  t  r  r 

£q.  (2-22)  for  the  mean  field.  Substituting  these  in  Eq.  ( 2- 
43)  a  d  simplifying,  w«  obtain 
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with  < * 


The  quantity  S  (k  )  is  the  transverse  newer  spectral  density 

C1 

at  the  interface.  The  expression  for  it  given  in  Eq.  (2  47) 
can  be  further  simplified  by  performing  the  integration.  We 
have 


i  >  Int<  >©  and  k  '  being  given  in  Eq.  (2-25) 
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Mow  applying  the  results  of  Appendix  8  relating  -  * ,  the 
backscattering  coefficient  directly  to  S  (k  ),  we  have 
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when  »k-sin-a  .  Noting  that  when  <  «k  we  have  <  **k  • 
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where  k’  <* 
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This  i.s  the  final  expression  for  the  backscattering  co¬ 
efficient.  Because  of  the  scalar  treatment  of  the  problem  the 
constant  y  can  not  be  evaluated.  The  situation  will  be  reme¬ 
died  in  the  next  section  where  the  same  techniques  os  used 
here  are  applied  directly  co  the  more  complex  vector  case. 
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III.  MULTIPLE  SCATTERING  OF  A  HORIZONTALLY 
POLARIZED  WAVE  BY  A.  RALF-SPACE  OF 
UNIFORMLY  DISTRIBUTED  DISCRETF:  SCATTERERS 


The  problem  of  multiple  scattering  cf  a  horizontally 
polarized  plane  wave  is  considered  in  this  chapter.  The  wave 
is  assumed  to  have  unit  amplitude  and  a  harmonic  time  depen¬ 
dence  of  exp(jwt).  It  is  incident  at  angle  •?  on  the  half¬ 
space  occupied  by  the  scatterers.  We  model  the  scatterers  by 
dielectric  spheres.  As  in  Chapter  Two  we  assume  that  the 
spheres  are  identical  and  randomly  distributed.  The  size  of 
the  spheres  is  again  taken  much  smaller  than  the  wavelength 
of  the  incident  wave  (Rayleigh  scattering) . 

The  main  quantity  of  interest  in  this  chapter,  as  in 
Chapter  Two,  is  the  backscatter ing  coefficient,  In  this 

chapter,  in  contrast  with  Chapter  Two,  the  scattering  coeffi¬ 
cient  la  evaluated  directly  in  terras  of  the  geometric  and 
electrical  properties  of  the  spheres.  This  allows  us  to  find 
an  explicit  expression  for  :*  in  terras  of  crown  quantities. 

In  determining  the  expression  for  •  *  we  follow  basically 
t he  s am*  steps  at  in  Chapter  Two . 

3-1  PROBLEM  FORMULA! ION  AND  FUNDAMENTAL  EQUATION 
A  horizontally  polarized  plan#  wave  is  incident  upon  the 
half  space*  at  random  spheres  at  an  angle  ? .  Tho  incident 
e l ec trie  field  £ n ( r ,  la  g i ve n  as 


x  s 


-  1  9  - 


*  e 


■  3  - 1 ) 


The  incident  field  is  scattered  by  the  dielectric  spheres,  and  thus 


gives  rise  to  a  scattered  field,  E^'r).  The  geometry  of  the 
problem  is  basically  the  same  as  the  scalar  problem  shown  in 
Fig.  2-1. 

We  now  derive  the  fundamental  equations  for  the  multiple 
scattering  of  an  electromagnetic  wave  in  the  Rayleigh  limit. 

As  ve  have  shown  in  Appendix  A,  the  spheres  can  be  represent¬ 
ed  by  an  equivalent  dipole  current.  From  electromagnetic 

theory  we  know  that  a  current  distribution  J(r')  produces  an 

( 9 ) 

electric  field  given  by 

Sfl  (r)  * ^  Gq  ( r-r ' )  .  Jg  (  r  '  )  dr  ’  (3-2) 

V 

where  G^ir-r')  is  the  free  space  dyadic  Green's  function,  and 

J- (r‘)  the  total  current  distribution  in  volume  V. 

-0 

The  free  space  dyadic  Green's  function  satisfies  the  fol¬ 
lowing  equations 

(  3  -  3  ) 

(  3-4) 

where  I-a  a  a  a  is  the  unit  dyadic  and  g,,  t  r)  is  the 
-  -x-jc  -y-y  —2  --  u  0  ••• 

free  space  Green's  given  by  Eq .  {2 -2a)  . 

The  total  current  J\  r’  -  is  given  by 
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Eq* .  (3-9)  and  !3~10>  are  the  fundamental  equations  of 
multiple  scattering  of  An  ■.iectrcmagnei.ic  wave  from  a  random 
collection  of  dielectric  spheres  vhich  are  small  compared  to 
wavelength . 


3-2  LV3lfWTIOH  If  THE  MEAN'  FIELD  OF  E(r) 

War  now  proceed  to  find  an  approximate  equation  for 
<E(r)>.  Following  the  same  steps  as  in  the  derivation  of  Eq. 
(2-19)  in  Chao ter  Two  we  have 


<E  ( r )  -  Eo(r)+:c/dr  5Jr-r  ).<E!r  )>  (3-11) 
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where  c  is  the  density  of  the  scatterers. 

Operating  cn  both  sides  of  Eq .  (3-11)  by  "x~x  s>  have 
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Note  that  the  incident  wave  i r !  satisfies  the  wave  equation 


m  free  space 


( . x7 x-x 
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9y  using  Eqs .  (3-13)  And  ; 3-35,  Eq.  (3-12)  becomes 
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By  using  Eq.  (3-11)  we  have 
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hMc  car.  writ®  this  43 


7x?x<E(r) > 
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(3-16) 


whore 
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W®  see  that  <E  { r ) >  satisfies  the  wave  equation  in  a 
■"continuous  tod .  uau“  in  which  the  propagation  constant  dopands 
upon  the  scattering  coefficients  and  density  of  the  scatter- 
er*.  Thu*  the  problem  of  finding  the  average  value  of  the 
electric  field  has  been  essentially  reduced  to  finding  the 
transmit  red  f-.eld,  in  a  half-space  (t>0)  of  "continuous  isedi'affi 
Carrying  out  this  calculation,  we  find  that  the  Man  wave 
in  the  equivalent  raediun  is 

1 (k  x~k  *2) 

*  E  ( r 5  »  -  T«’  *  2  a  ,  z>0  (3-18) 


where 
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and  k  and  k  ars  given  by  £q,  (2-24)  ar»d  (2-26)  respecti vely . 
For  partis® ter  a  values  typical  of  vegetation,  we  will  find  that 
:J<<1  which  implies  the  T*l. 

3-4  EVALUATION  OF  THE  TRANSVERSE  SPECTRAL  .DENSITY  AND 
BACK SCATTERING  COEFFICIENT 

In  this  section  we  evaluate  the  transverse  power  spectral 
density.  The  procedure  differs  slightly  from  the  previous 
section.  Here  the  Fourier  transform  of  the  scattered  field  is  carpeted 
first  and  then  the  spectral  density  is  obtained.  Previously 
we  computed  the  correlation  of  the  scattered  field  and  then 
tr  ana  t o  rrmd  i  t . 

The  single  scattered  field  in  the  equivalent  medium  is 
obtained  by  modifir.g  £q.  (3-2).  We  have 
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where  G  and  J  have  replaced  Gg  and  J4g .  Here  G  is  the  dyadic 
Green's  function  in  the  equivalent  .medium.  It  is  given  by 
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with  k  _  being  given  in  Eg.  { 3-17) .  The  induced  cur ren 
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Next  w®  take  the  trt  saverao  Fourier  transform  of  the 
scattered  field  {the  transform  was  defined  in  Eq .  (2-44)!. 
Noting  the  fact  that  Eq.  (3-21)  is  a  convolution,  we  have 
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where  E  fk  ,s)  x»  the  Fourier  transform  of  E  (r) 
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If  wf  now  integrate  the  second  term  in  Eq,  (3-28)  by  part® 
to  raraove  ti*e  derivative*  with  respect  to  x  and  y„  we  obtain 
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The  Fourier  transform  of  the  scattered  f^eld  can  novr  b® 
computed  explicitly.  By  using  Eqs.  (3-26)  and  (3-31)  in  Eq, 
(3-25),  we  obta  •it  *  i 
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where 
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tfext  w#  resnov®  the  apectdlar  field  from  consideration  by  de¬ 
fining  the  fluctuating  f ield  as  follows: 
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Now  %m  forns  the  dot  product  of  Ej  with  its  conjugate  and 
average.  Letting  V,  N-~»  such  that  N/V»p,  we  find 


(3-37) 


Substituting  Eq.  {3-35}  into  Eq.  (3-37) ,  it  follows  that 
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where  the  mean  field  from  Eq.  (3-18)  has  been  used. 

Since  we  require  the  spectral  density  in  the  z»0  plan® 
we  first  use  the  fact  that 
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in  Eq.  (3-38)  anti  then  set  z,*z,-0.  We  obtain 
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front  appendix  S,  we  have 
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Since  there  is  no  k  component  of  the  incident  waive,  we  have 
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Now  putting  Eq.  <  j: — 4 0 )  into  Eq .  (3-41),  we  obtain 
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By  using  Eq.  (3-29)  in  Eq.  (3-43)  and  evaluating  the  integral. 
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with 
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In  the  study  of  vegetation  such  as  forest,  canopies,  w# 
:ind  that  oa<<l-  Then  Eq.  (3-45)  can  be  written  approximately 
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Using  this  in  Eq.  (3-19),  we  see  T  si,  i.e.  ,  there  is  very 
little  reflection  at  the  interface.  Eap.loyi.no  this  result 


and  the  small  pS  approximation  given  in  Eg .  (3-46),  the  for¬ 
mula  for  aa  given  in  Eq.  (3~44)  can  b«  substantially  simpli¬ 
fied.  Me  find 
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- TTr^l - 
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This  is  our  final  result.  We  see  that  it  is  independent 
of  the  density  p  to  first  order  in  p&.  The  angular  variation 
is  a  cosine.  This  angular  dependence  corresponds  to  the  third 
empirical  model  proposed  by  Clapp  ,  however,  Clapp  did 

not  determine  the  multiplicative  constant  as  we  hat's  for  the 
case  of  small  spheres. 

Various  limitations  on  the  above  formula  should  be  point¬ 
ed  out:  first,  the  Foldy  closure  assumption  which  allowed  us 
to  obtain  an  equation  for  the  mean  «f<*ve  is  moot  likely  only 
good  when  8 c <  <■  1  ?  second,  the  backacatter  angle  must  be  bounded 
away  from  grazing  or  the  approximation  in  Eq.  (3-46)  will  not 
be  valid;  and  third,  the  spheres  must  have  sufficient  lo»s  so 
that  the  mean  wave  does  not  penetrate  toe  far  into  the  medium. 
This  insures  the  appl iabil ity  of  the  Born  approximation  to 
compute  the  backscatter. 
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IV.  EMPIRICAL  TECHNIQUES  AND  NUMERICAL  RESULTS 

In  this  chapter  we  compare  our  results  to  an  empirical 

( 12 ) 

model  developed  by  Attema  and  Ulaby  .  The  comparison  aids 
in  the  physical  interpretation  of  the  results  obtained  in 
Chapter  III.  Following  this,  we  evaluate  the  backscattering 
cross-section  obtained  in  Chapter  SIX  for  various  parameters. 
Plots  of  o’*  versus  the  angle  of  incidence,  3,  are  presented 
for  different  value*  of  frequency  and  sphere  radii  a. 

4-1  EMPIRICAL  MODELING 

We  have  studied  the  scattering  properties  of  vegetated 
terrain  by  treating  the  target  as  a  collection  of  lossy  die¬ 
lectric  spheres  and  deriving  the  backscattering  cross-section 
directly  from  Maxwell's  equations  and  the  statistics  of  the 
medium.  Attema  and  Ulaby  treated  the  same  problem  by  a  usore 
empirical  approach  which  we  will  describe  below. 

They  assumed  the  vegetation  could  be  modeled  by  a  layer 
of  water  droplets  having  thickness  h.  The  basic  geometry  of 
the  model  is  shown  in  Fig.  4-1.  To  keep  the  model  as  simple 
as  possible,  the  following  assumptions  were  made:  first,  the 
water  droplet  cloud  representing  the  vegetation  consisted  of 
identical  water  particles,  uniformly  distributed  throughout 
the  lay* r;  second,  only  “single  scattering"  was  considered. 
Her®  single  scattering  was  taken  to  mean  single  scattering  by 
the  effective  or  average  wave  in  the  medium. 

The  reflectivity  factor  or  radar  cross-section  per  unit 
volume  and  the  power  attenuation  coefficient  per  unit  length 


Fig.  4-1  B®an  Incident  Upon  Slab  of  Vegetation 

were  expressed  as 

--  •»  c  ( 4  - 1 ) 
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where  t  is  the  number  of  water  particles  per  unit  volume,  j 
is  the  radar  back, scattering  cross-sect ion  of  &  single  droplet 
and  0  i*  tii*  total  attenuation  crons  -section  for  a  single  par 
tide . 

To  calculate  the  average  backscattered  power  due  to  the 
incident  radar  bean,  the  eonfcr ibut ions  by  all  particles  in  th 
beam  are  sunaned,  taking  into  account  th*  two  -wav  attenuation 
bv  the  vegetation  .Laver  between  the  scattering  pa  rticle  end 
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the  vegetation  surface.  Referring  to  Figure  4-1,  the  inci¬ 
dent  w* ve  is  assumed  to  be  a  plane  wave  with  power  density  S 
confined  to  an  ideal  cylindrical  beam  illuminating  the  surface 
at  an  angle  o£  incidence  9.  Let  the  illuminated  area, 
be  defined  as  the  intersection  of  the  beam  with  the  horizon¬ 
tal  plane  at  the  top  of  the  vegetation  layer.  The  incident 
power  will  then  be  given  by: 


SAj 11cose 


(4-3) 


and  the  average  back  scattered  pc  we ::  <«>r>  is  then  found  to  be 
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Consequently, 
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This  is  the  final  formula  they  used  to  relate  to  the  me¬ 
dium  cross  .section#  and  density. 

For  our  purpose#  the  formula  can  be  simplified  by  first 
assuming  uh  '  >  1  ?  thus  Eq .  (4-5)  becomes 


Mow  we  specialize  cons iderations  to  Tlayleigh  scattering,  i.e.  , 
we  assume  the  wavelength  is  large  compared  to  the  droplet  size . 
If  w*  further  assume  the  droplets  are  spheres  of  radius  a , 
then  the  backseat taring  cross -sect ion  of  an  individual  sphere 


is  given  approximately  by 
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where  K  has  been  defined  in  Appendix  B.  Next,  we  write  Q  as 
a  sum  of  an  absorption  cross  sectionj  Q  ,  and  scattering 

A 

cross-section,  Q 


Q  -  0a  *  Qs 


(4-8) 


In  the  Rayleigh  limit  for  spherical  scatterers,  we  have 
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For  the  X  band  region  of  the  spectrum  under  consideration,  the 

attenuation  cross-section  is  much  larger  than,  the  scattering 

cross-section;  i  . e . ,  Q  -  • Q  .  Thus 

s  a 
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Now  using  L'qs,  (4-2),  (4-7)  and  '4  —  11}  in  Eq .  ( 4-6)  ,  we  ob- 
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This  is  exactly  the  expression  obta ined  for  -.?*  in  Chapter  III 
Thus  we  see  in  the  Rayleigh  limit  for  spheres  w«  can  be  *a- 


tisffied  that  the  empirical  technique  and  the  analytic  vector 
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method  cive  the  same  results. 

4-2  NUMERICAL  EVALUATION  OF  c° 

In  this  section  the  backsc&ttering  cross-section  as 
given  in  Eq.  (4-125  is  evaluated  fcr  several  sphere  sizes 
and  frequencies  as  a  function  of  incidence  angle,  3.  We  as¬ 
sume  that  the  spheres  are  droplets  of  water.  The  complex 

(16) 

permitivity  of  water  as  given  in  Peake  is 


£ 
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where  \  is  the  wavelength  given  in  centimeters. 

The  plots  are  shown  in  Figs.  4-2  and  4-3.  The  frequen¬ 
cies  used  were  9  and  15  GHz  respectively.  An  examination  of 
the  plots  shows  that  the  general  behavior  of  the  plots  la  in 
agreement  with  experimental  data1  if  spheres  in  the  order 
of  a  few  millimeters  are  chosen.  To  obtain  more  definitive 
information  the  model  in  Chapter  III  would  have  to  be  gener¬ 
alised  to  consider  discs  and  cylinders.  This  work  is  in 
progress . 


Backscatter inrj  Coefficient  versus  Angle  of  Incident  at  15  Gtiz 
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V.  DISCRETE  ONE  DIMENSIONAL  RANDOM  SLAB 


The  proceeding  two  chapter*  have  been  concerned  with 
beckscattering  from  three  d ut«*n*ion*i  scatterers  that  were 
randomly  distributed .  The  calculation  of  the  faackscattering 
coefficient  was  perfc  -aed  by  employing  Foldy's  approximate 
method  to  calculate  the  asar,  field.  In  this  chapter  we  will 
compare  the  Foldy  approximation  for  the  average  dielectric 
constant  to  the  expression  obi  ,n*d  in  the  thermodynamic  lim¬ 
it,  that  is,  in  the  limit  when  particle  size  becomes  small 
and  particle  density  large.  To  accomplish  this,  we  will  re¬ 
strict  our  attention  to  one  dimensional  discrete  model*. 

Here  the  medium  is  composed  of  identical  dielectric  slab* 
which  are  distributed  probaIi*ticaiiy .  The  results  obtained 
should  carry  over  to  the  three  dimension  case,  however,  this 
work  has  not  been  completed  at  th  s  time . 

The  saethods  used  in  this  chapter  differ  substantially 
from  those  employed  previously.  First , identical  slab* 
of  width  "a*  are  chosen  so  that  they  obey  a  renewal  process 
which  will  be  described  in  »or*  detail  in  the  next  section. 


This  procss 
9  lab*  never 
distributed 


he*  the  important  property  that  overlapping  of 
c cur*  If  on*  choose*  the  slab  position  to  to* 

n  a  Poisson  manner ,  then  overlapping  will  occur. 
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Using  a  process  without  overlapping  is  desirable  sine©  this 
correrponds  to  more  physical  situations  of'  int#r®at. 

Following  tit®  choice  of  the  renewal  process, 
differential  equations  are  found  for  the  probability  of  occur¬ 
ence  and  nonoccurence  of  a  slab.  The  solutions  to  these 
equations  enable  or©  to  calculate  the  average  dielectric 
constant  of  the  random  media's  exactly. 

In  the  remainder  of  the  chapter,  the  renewal  process  is 
used  to  drive  a  nonlinear  system  of  ordinary  differential 
equations  having  solution  u(tj.  A  diffusion  equation  is  de¬ 
rived  for  the  probability  density  of  u(t)  and  the  solution  to 
this  equation  is  obtained  in  the  thermodynamic  limit.  Final¬ 
ly,  these  general  results  are  applied  to  the  one  dimensional 
slab  and  compared  to  Foldy’s  results. 


5-1  SLAB  RENEWAL  PROCESS 

In  thi a  section  we  will  describe  the  random  process  that 
we  wi.il  use  to  mode  1  the  one  dimensional  random  medium.  The 
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where  r.  i.a  the  free  space  permittivity  and  is  a  constant 
da  termini  ng  the  strength  of  the  fluctuations.  The  process 
tiz)  takes  on  values  0  or  1  if  z  is  outside  a  slab  or  in¬ 
side  a  slab  respectively. 

The  probabilistic  character  of  the  process  r (*)  will 
now  be  specified.  Associate  the  random  variable  Z ^  with  th* 

fJl 

trailing  edge  of  the  t  slab.  Rather  than  specify  th©  Z 
directly,  we  form  the  difference 
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where  2  *  0 .  The  W  ar®  vast  the  distance  between  trailing 

o  1 


edges  of  adjacent  slabs.  Me  assume  that  the  W  arc  independent 


hi a  definition  of 
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identically  distributed  random  variables, 
the  W(  forces  the  process  r lz>  to  o«  a  renewal  process 
choose  the  distribution  of  the  w.  to  be 
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where  u • w •  is  the  step  function ,  for  small  "a " ,  the  parameter 
■  can  be  interpreted  as  the  number  of  pulses  per  unit  length 

or  the  puls®  density . 

This  assumed  distribution  has  two  important  properties . 
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“lapped  the  dielectric  constant  would  double  in  the 
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second  property  of  Eq.(5-*3),  is  that  for  small  "a”  the  dis¬ 


tance  between  slabs  is  approximaualy  exponentially  distributed 


This  implies  that  the  slabs  ar®  randomly  distributed  which  is 


physically  appeal  m  9 ,.  It  the  2^  had  been  chosen  to  be  uni¬ 


formly  distributed  at  the  outset,  then  the  process  would  have 


been  Poisson.  This  is  &  usual  assumption  made,  however,  it 


leads  to  an  overlap  problem. 


Since  the  process  r(2)  has  now  beer,  defined,  some  of  its 


properties  can  be  calculated  from  its  probability  density  ?  (z) 


where 


p_  ( z) 
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Here  p  (2  is  the  probability  that  the  ccint  z  is  outside  a 

o 


slab  while  p  1  (  z )  is  the  probability  that  z  is  located  inside 


a  slab.  We  see  that 
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Thus  *  knowledge  of  p. (z)  will  allow  us  to  calculate  the  average 


dielectric  constant.  Using  the  method  of  conservation  of  pro- 
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babilitv  one  can  show  that  p^ii!  and  p^iz}  satisfy  the 
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Eq  •  (5-7)  can  be  solved  by  employing  Laplace  transform  tech¬ 
niques,  Denoting  5  (s)  as  the  Laplace  transform  of  p  (z)  , 

o  o 
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the  solution  to  Eq.  (7)  becomes 
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Now  employing  the  inverse  Laplace  transform  p  (z)  becomes 
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Now  using  Eqs.iS-5)  (5-65  ,we  find 
«■  c  (  z )  >  * 


(5-12) 


In  general  the  mean  will  be  a  function  of  z  and  thus  the  pro¬ 
cess  is  net  stationary;  however,  an  examination  of  the  mean 
for  large  z  shows  it  approaches  a  constant  value.  By  using 
the  final  value  theorem,  we  find 

1  im  p  *  z )  <«  i  .*.  ,t.  s  ;  i  s  5  *•  .-it. — -  (5-13) 
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A  n*or<»  detailed  calculation  snows  uiat  tan  constant  mean  is 
approached  w.thm  a  few  pulse  distances .  Thus  tne  nonstaci.cn- 
arv  character  of  the  mean  is  ;  ust  due  to  the  boundary  effects 
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at  i*o. 


Before  concluding  this  section,  we  will  discuss  the 
thermodynamic  limit.  In  this  limit  we  let  the  puls*  width 
a~o  as  the  pulse  density  V— .  We  do  this  in  such  a  way  that 
Aa»constant*i.  In  this  case,  the  medium  approaches  a  non- 
trival  limit  which  we  call  the  macroscopic  ot  thermodynamic 
limit.  Performing  this  limiting  operation  on  Eq.  (5-12) ,  we 
find 
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°i 

Thus  we  see  the  mean  has  a  constant  value  for  all  z  greater 
than  zero.  w«  shall  return  to  this  behavior  at  a  later  point. 

5-2  NONLINEAR  SYSTEM  DRIVEN  BY  RENEWAL  PROCESS 

We  will  now  divert  our  attention  from  the  specific  study 
of  the  one  dimensional  problem  and  obtain  a  result  that  will 
be  needed  in  the  next  section.  Consider  and  Nth  order  system 
of  nonlinear  ordinary  differential  equations  driven  by  a  re¬ 
newal  process  r ! z )  i.e. 
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where  j(u-u  )  is  an  N  dimensional  delta  function.  Here 
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u(Cju,z)  is  the  solution  to  the  final  value  problem 
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and  ■'(C,u,2)is  the  Jocabian  of  the  transformation  from  u  to 
u ,  i  .  c  . 

30. 


J  ( *; ,  a,  z)  *  det 


'5-191 


Although  the  aquation  is  complicated,  one  must  realize  it  is 
an  exact  equation  with  no  approximations  for  p  (u,:i  and  the 
equation  is  deterministic. 

Again  employing  conservation  of  probability  arguments , 
an  additional  equation  can  be  obtained  that  relates  p, (u,2) 
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where  p,  i s  defined  in  F.o.  {5-1.65  and  (z)  .  is  z  if  z>o  and 
zero  if  x^p. 

We  would  now  like  to  see  how  these  equations  reduce  in 
the  thermodynamic  limit,  i.e.,  as  a-o  and  v~«  such  that  a>*a 
Expanding  Eqs .  (5-17)  and  (5-20)  in  a  power  aeries  in  “'a",  re¬ 
placing  \  by  s/a  and  keeping  only  dominant  terms,  w®  obtain 
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The  probability  density  with  respect  to  just 
by  summing  out  the  r  variable.  We  have 
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can  be  obtained 
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»y  using  Eqs  .  (5-21)  and  (5-22)  an  equation  for  p  <  u ,  x ) 
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If  we  now  make  the  important  assumption  that  £  depends 
linearly  upon  the  renewal  process  r  then  we  car.  write 
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this  result  in  Eq .  (24) 
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3y  direct  substitution,  it  can  be  shown  that  the  solution  to 
Eq.  (5-26)  is 
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Thus  w«  see  the  probabi 
is  a  delta  function,  i. 
value  u(2i.  Therefore 
that  this  deterministic 

a  s  t  h «(  r  a  ndons  solution. 
Thus  in  this  iin.it  the 
There  are  so  many  parti 
sees  the  averse,®  effect 


lity  density  in  the  thermodynamic  limit 
a.,  at  each  point  z,u  only  takes  on  one 
the  solution  is  deterministic .  Mote 
solution,  u( 2) ,  obeys  the  same  equation 
Uvz)  ,  Eq.  (5-L5)  with  r  replaced  by  <  r>  . 
process  acta  in  a.  deterministic  way. 
cles  per  unit  length  that  the  system  onl 
of  the  particles .  This  result  will  be 
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random  medium.  We  shall  calculate  the  average  re¬ 
flected  wave  and  compare  these  results  with  those  obtained 
by  the  Fcldy ‘ 3  technique. 

Consider  a  scalar  problem  having  field  component  u  which 
obeys  the  reduced  wave  equation: 


dz 


k  <s  (  z )  u 
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—  «•  <  2  <  ** 


(S-29: 


when  k  is  the  free  space  wave  number  of  the  medium  and  s(z) 
is  the  dielectric  permittivity.  The  domain  of  the  variable 
permittivity  is  confined  to  o<z<L.  Outside  this  region  c ( z > 
•takes  on  a  constant  value  one.  Inside  the  slab  region  the 
dielectric  fluctuations  ere  mooelled  by  a  renewal  process  as 
follows : 

1  i.  ,  z>i, 

(5  -  JO 


i  ( z) 


2  '! 


*>L 

O  <  Z  <  L 
7.  <0 


where  r  {?.)  has  been  defined  earlier.  We  assume  a  plane  wave 
is  normally  incident,  upon  the  slab  from  the  right  and  a  re¬ 
flects  and  transmitted  wave  is  generated  in  the  homogeneous 


region  to  the  right  and  left 


:he  slab  respectively.  This 


.3  shown  i  n  rid.  5  - 1 ... 

Tne  field  m  the  homogeneous  regions  can  be  written  as 
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where  "  L>  is  the  reflection  coefficient  of  the  slab  and  T.L) 

1 3  the  transmission  coefficient  .  w»>  have  exol  icitl v  indicated 


that  both  r  and  T  are  functions  of  the  slab  thickness. 

Since  we  want  to  draw  upon  the  results  of  the  last  section, 
the  boundary  value  problem  of  £q.(5-29)  will  have  to  be  con¬ 
verted  to  an  initial  value  problem.  This  can  be  accomplished 
by  looking  at  the  reflection  coefficient  "(L)  rather  than  the 
field  u(z) .  One  can  derive  a  Rxcatti  equation  for  this  re¬ 
flection  coefficient.  It  is 
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We  see  it  is  a 

first 

order  nonlinear 

ordinary 

dif  ferential 

equation.  The  problem  is  now  an  initial  value  problem  since 
the  condition  in  Eq .  (5-33)1.3  only  given  at  one  point. 

The  problem  has  almost  been  cast  in  the  form  of  the  last  sec¬ 
tion  except  for  the  fact  that  Eq .  (3-32)  is  complex .  To  trans¬ 
form  Eq .  (5-30)  into  a  real  form,  we  let 
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where  is  the  amplitude  of  7  and  j  is  its  phase.  Plugging 
this  into  Eq .  and  equating  real  and  imaginary  parts,  we 

find 
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replaced  by  L  and  with 


In  addition,  since  the  system  is  linear  in  r(L),  the  results 
derived  for  the  thermodynamic  limit  hold  for  Eq.(5-35).  Bas¬ 
ically  this  says  that  when  a  -o  and  v *»  such  that  al«a  the 
solution  for  r  and  i  is  deterministic  and  can  be  obtained  by 
replacing  r(L)  by  <r(L)>.  Physically  this  says  that  in  the 
thermodynamic  limit  the  solution  to  the  problem  can  be  ob¬ 
tained  by  replacing  the  statistical  medium  rtz!  by  its  average 
•r(z)  Thus  the  medium,  can  be  replaced  by  one  having  an 
average  dielectric  constant:  of 
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An  analysis  of  the  one  dimensional  problem  by  the  1'oldy 
technique  leads  to  the  following  result,  for  the  average  di- 
e  1  e  e  t  r  i  co ns  t  a  n  r.  : 

«  z  '■  >  .  i  ■.  1  ♦  ■  i  i  (  S-  3  7  ) 

i  o 

Thus  we  see  that  the  result  s  from  the  Fold’/  technique  over¬ 
lap  with  t  he  thermodynamic  1 imit  results  when  * •  -  1 .  It  also 
indicates  that  the  Poldy  method  is  not  correct  when  •  «•* ,  a  *o 
aid  • a  is  moderate  *  r  large. 

Although  the  analysis  at  present  has  net  resolved  the 
nature  of  the  di f  fe rents  between  the  "anal  tie  the iry“  and  the 
i  unspo  t.  t  heory "  ,  the  existence  n*  an  exact  equation  for  t  h  e 
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probability  density  should  prove  to  ha  a  useful  tool  in 
the  future- 


XV.  CONCLUSION  AND  RECOMMENDATIONS 


Wo  have  developed  and  analysed  a  model  for  a  half  spaca 
of  vegetation.  This  model  views  th®  vegetation  as  a  collection 
of  lossy  dielectric  spheres*.  Within  the  constraints  of  this 
model  the  hackscattering  cross  section  ha®  been  related  to 
c  here  sire  and  dielectric  constant.. 

In  Section  II  the  scalar  problem  was  analyzed  by  the  Foldy 
technique.  This  section  served  the  purpose  of  illustrating 
the  technique,  however,  because  of  the  scalar  nature  of  the 
problem  a  complete  relationship  between  th®  .medium  parameters 
and  hackscattering  cross  section  could  not  oe  obtained.  In 
Section  III  this  defect  was  remedied  by  applying  the  Foldy 
technique  to  the  complete  electromagnetic  problem.  The  case 
of  horizontally  polarized  waves  was  treated  and  an  expression 
was  obtaxned  relating  the  radar  crocs  section  to  the  m<wi i urn 
parameters.  ±his  expression  was  found  in  agreement  with  the 
angular  dap tndence  of  Clapp's  third  model.  An  interesting 
feature  of  this  expression  is  that  it  was  independent  of  den¬ 
sity. 

In  Section  IV  the  results  of  Section  ITT  were  compared 
with  the  data  presented  by  Atteraa  and  Uiaby,  The  agreement 
was  found  to  be  quit®  good  for  water  droplets  of  1  to  5  rwn  in 
the  X  band  region  of  the  spectrum.  In  addition  a  connection 
was  established  between  the  more  empirical  aralysir  of  Attests, 
and  Uiaby  and  the  systematic  procedure  employed  in  Section 
Itl„  The  connection  showed  that  the  method  we  employed  ha#, 
it,  effect,  assumed  that  the  scattering  cross  section  of  an 


individual  particle  is  small  compared  to  the  absorption  cross 
section.  The  validity  of  this  assumption  in  the  frequency 
range  of  interest  i.®  shown  in  Figure  4-4. 

Finally  in  Section  V,  we  explored  some  fundamental  limi¬ 
tations  cf  the  Foidy  technique.  There  wo  saw  by  employing  a 
one  dimensional  mode,  that  the  Foidy  approximation  only  gives 
correct  results  when  the  product  of  slab  width  ar.d  density  are 
smal 1 . 

In  view  of  the  success  of  the  method  we  used  in  relating 
physical  parameters  of  the  ctcdi.-m  re-  har.kacattor.lng  croeu  sec¬ 
tion,  we  make  the  following  rt  courrendatif,  u . 

1.  Analyse  a  medium  represented  by  either  having  a  ran- 
doss  local  lex*  arid  a  probabilistic  angular  orient  At  Ion. 
Sine®  the  disc®  would  represent  leaves,  the  wave¬ 
length  would  b«  of  the  order  of  the  ace  turn.  A  no  object 
Twersky’s  technique  (a  general irutd.oo.  of  Foidy**  me¬ 
thod  which  remove®  the  Rayleigh  assumption)  would  be 

employed.  The  results  would  be  compared  to  'those  of 
(16) 

Dn  who  calculated  hackscattering  from  leaves  by 

a  different,  toot «  heuristic,  technique, 

2,  By  using  the  technique  employed  in  recommend at ion  1., 
analyte  the  effect  of  ground  lying  unde."  the  veoeta- 
tier.  This  requires  the  analysis  of  the  slab  problem 

3  investigate  the  effects  of  depolarisation  ,  S*a  how 
they  compare  to  tha  results  of  continuous  random  me¬ 
dium  nv’JcUsl*. 
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Ap  endix  A 

Induced  Current  on  a  Dielectric  Sphere 

When  the  site  of  a  dielectric  sphere  is  much  smaller  than 
the  wavelength  of  the  radiation  incident  upon  it,  the  electric 
field  within  and  near  the  sphere  behaves  like  a  constant  elec¬ 
trostatic  field.  From  elementary  electrostatics,  we  know  the 
scattered,  field  can  be  represented  by  a  dipole  located  at  th® 
center  of  the  sphere.  In  this  appendix  we  will  calculate  the 
equivalent  current  J  (r)  .  This  will  be  derived  from  the  in- 
duced  dipole  moment  as  given  by  Stratton. 

Let  us  assume  there  is  an  electric  field  £(r)ay  incident 
upon  a  sphere  of  radius  a  whose  center  is  located  at  r»r ' . 

If  the  sphere's  radius  is  small  compared  to  wavelength  then 

(9) 

the  induced  dipole  moment  density  P  (r )  is  given  by 

?(r)  *  P(r )  ay  »  <’TtQka3E  £r ' )  5  (r-r ' )  fly  (A-l) 

whore  K« ( £  -]»/{£  *2)  and  e  is  the  relative  dielectric  con- 

XT  XT  X 

stant  cf  the  dielectric  sphere. 

The  equivalent  charge  density  o  (r)  used  to  represent 

the  scattered  field  is  then 

;>  (r)  -  - 7 . ? ( r)  (A-2) 

#q  -  - 

*  -4rr  t3ka'E<  r  ' )  '*  (  x-x  ’  )  5 '  (y-y  *  !  5  (  z-z  • }  (A- 3} 

when  5 1  (z)  is  the  derivative  c-f  Mr;  with  respect  to  z.  From 
the  continuity  equation 

-54- 


e-sis 


( A-4 ) 
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(A-6) 


Thus  we  can  replace  each  sphera  whoa©  radiua  is  amall  compared 
to  wavelength  by  equivalent  dipole  current  as  given  in  Eqa. 
(A-5)  and  (A-6) . 
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Appendix  B 


Relationship  Between  Backscattering  Coefficient  and  Transverse 

Spectral  Density 


In  both  Chapters  II  and  III  w®  have  used  the  relationship 
between  transverse  power  spectral  density  5 { k ^ )  and  o®.  In 
this  appendix  we  will  deriv®  that  relationship. 

Assume  that  a  vector  component  of  the  scattered  field, 
say  ( rt , 0 )  is  known  on  the  boundary  (r-0  plane)  between  the 
random  medium  (z>0)  and  free  space  (z<0)  .  Then  the  far  field 

in  the  free  space  region  due  to  a  region  A  on  the  plane  xm0 

1 1 7  \ 

is  given  by  the  Kirchhoff  diffraction  formula  .  It  is 

"jk0r 

jk.cosSe  f  jk  x'sine 

♦.<r) - !Vr - /v-:;-0'*  <»-i» 


where  9  is  measured  with  respect  to  the  normal  to  the  £™0 
plana  as  is  shown  in  Fig.  2-1  and  r*  j  rj  » rxT+yr+jt'* . 

The  backscattering  coefficient  is  defined  as 
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where  If  is  the  average  intensity  at  the  receiver  of  ,mo>  fluc¬ 
tuating  scattered  field,  i.e.,  X  •*<$* (r)  $5(r)  >  wnere  $  -<$  >, 

and  is  the  field  intensity  incident  upon  the  illumina¬ 
ted  area  A.  Assuming  an  incident  field  of  the  form 

jk^ <  Kains-zooefi ) 

>  j  (£)■  w  e  '  (B~3) 
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then  we  have  I ^*! (r) ! **1. 

Th«  average  at  th«s  receiver  from  Eg.  (B-l)  is  given  by 
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tAeare  ( r " ,  0)  *$  ( r ' ,  0)  -<0- (rj,  #  0)  > .  Using  Eqs ,  { a~*  l )  and  (8-4) 
r  v  •  —t  r 


in  Eq.  (B-2),  we  have 
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W«  will  now  use  the  transvei  ?e  Ft  ur  ier  transform  of  ?f  (r^ ,  0) 
and  .its  conjugate.  Denotin  the  'm  "ora  of  ,  .,„fr^./0)  by 
(k t, 0? ,  we  have 
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Us in  g  th is  in  Eq .  (B-5. , 
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whet »  S(k^)  is  the  transverse  power  spectral  density  of  the 
interface  field.  Putting  Eq.  (S-8)  into  Eq.  (B-7)  and  making 
the  change  of  variables 
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we  find 
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Illuminating  the  whole  plane  ( A •*« ) ,  the  fa i&l  result  1m 
obtained.  It  is 


k^coa*  9 

a •  «  — — k — —  S  C  k  ) 

4>rr  ~t 


(B-1.D 


where  k  « k ,, 
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